AO-A069  336 


UNCLASSIFIED 


COLORADO  UNI V BOULDER  DEPT  OF  MATHEMATICS  «r/6  15/7 

ON  QUALITATIVE  LANCASTRIAN  MODELS  OF  COMBAT  INCORPORATING  LOGI — ETC(U) 
MAT  79  J S MAYBEE  DAAG29-78-6-0113 


ARO-16105.1-A-M 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 


Unclassified ____ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWnn  Of  Enl><0 


REPORT  DOCUMENTATION  PAfct 

' i Ivt^ceIwIn 


4.  T\V^P^^ubtltto^^ 


^ READ  INSTRUCTION  S 
BEFORE  COMPLETING  FORM 


p-'V  m 


On  Qualitative  Lanchestrian  Models  of 
Combat  Incorporating  Logistics  . 


RECIPIENT'S  CATALOG  NUMBER 


OF  REPORT  A PERIOD  COVERED. 

Technical  r£ p' 


>-AUTMORf»?S=ii:= 

| q John  S.^Maybee  J 


PERFORMING  ORGANIZATION  NAME  AND  AODRESS 

University  of  Colorado 
Boulder,  Colorado 

(I.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

U.  S.  Army  Research  Office 
, P.  0.  Box  12211 

Research  Triangle  Park,  ..'C  27709 

14.  MONITORING  AGENCY  NAME  * AOORES57II  dlllmtont  from  Controlllnh  Office) 


— — — » — — . 

6.  PERFORMING  ORG.  REPORT  NUMBER 

S.  CONTRACT  OR  GRANT  NUMBER! •) 

Y{£  I /DAAG29-78-  G-jfi  1 ^ 

k!  

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  a WORK  UNIT  NUMBERS 


n 


May  79 

"number  cJT'p  AGES 

23 


1 


15  SECURITY  CLASS,  (ot  thlm  report) 

Unc lassified 


IS*  DECLASSIFICATION^  DOWNGRADING 
SCHEDULE 


DISTRIBUTION  STATEMENT  fof  thlm  Report) 

Approved  for  public  release;  distribution  unlimited. 


I 17.  DISTRIBUTION  STATEMENT  (of  tho  mb,  tract  mntmrmd  In  Bloch  iO.  II  dlllmrmnt  from  Report) 


D D C 

^ JUN  4 1879 


I ia.  supplementary  notes 

The  view,  opinion?,  and/or  findings  contained  in  this  report  are  those  of  the 
author ( s ) and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  documentation. 


[ 19.  KEY  WORDS  ( Continue  on  rovoroo  tide  it  nocootary  and  Idontity  by  block  number) 


Logistics 

Models  of  military  combat 
Lanchestrian  models 
Differential  equations 
Eigenvectors 


l4  y sj  o ^ 


u 


1 20 


A 


j - ■ - _ 

ABSTRACT  (Continue  on  roooroe  tide  It  nocoooary  and  Idontity  by  block  number) 


We  initiate  the  study  of  a family  of  models  of  military  combat  which  include 
classical  linear  Lanchestrian  models  as  a special  case.  These  models  introduce 
the  additional  concept  of  one  or  several  columns  providing  supplies  to  one  of 
the  forces.  It  is  assumed  that  the  other  force  is  defending  an  already  supplied 
position  or  is  prepared  only  for  a short  campaign.  The  models  are  designed  so  asj 
to  have  solutions  that  are  primarily  qualitative  in  character. 


tative  in  cnaracter.  j . 


oo M73 


EDITION  OF  I NOV  45  IS  OBSOLETE 


Unclassified 


litatlve  Lanchestrian  Models  of 


Combat  Incorporating  Logistics 


John  S . Maybee 
Department  of  Mathematics 
University  of  Colorado 
Boulder,  Colorado  80309 


This  work  was  supported  by  a grant  from  the  U.S.  Army  Research  Office 


Acoesai 


CRA&I 


Abstract 


We  Initiate  the  study  of  a family  of  models  of  military  combat  which 
include  classical  linear  Lanchestrian  models  as  a special  case.  These 
models  introduce  the  additional  concept  of  one  or  several  columns  pro- 
viding supplies  to  one  of  the  forces.  It  is  assumed  that  the  other  force 
is  defending  an  already  supplied  position  or  is  prepared  only  for  a short 
campaign.  The  models  are  designed  so  as  to  have  solutions  that  are  pri- 
marily qualitative  in  character.  They  have  the  form  of  a linear  system 
of  differential  equations  -r  = Au  + F(t)  where  u is  a vector  with 

> sr  sr  <fS, 

m X 3 components,  A is  an  m,x  m constant  matrix  and  F is  a given 
time  dependent  vector  with  m conponents.  The  differential  equations 
are  to  be  solved  subject  to  certain  stopping  rules.  When  we  say  that 
the  solutions  are  qualitative  in  character,  we  mean  that  the  mathematical 
form  of  the  solution  does  not  depend  upon  the  magnitudes  of  the  elements 
of  the  matrix  A but  only  upon  the  signs  of  these  elements.  We  intro- 
duce the  concepts  of  upper  and  lower  solutions  which  serve  as  bounds 
for  all  positive  solutions  of  the  problem.  We  obtain  detailed  results 
on  the  eigenvalues  and  eigenvectors  of  A and  explicit  representations 
of  the  upper  and  lower  solutions. 


1.  Formulation  of  the  Models 


In  this  paper  we  introduce  a class  of  Lanchestrian  models  of  com- 
bat Including  logistic  considerations  and  Initiate  the  detailed  analytic 
study  of  them.  It  turns  out  that  the  analysis  of  such  models  falls 
naturally  into  two  parts,  an  analytic  part  in  which  we  derive  explicit 
solutions  and  examine  their  fundamental  properties  and  a geometric  part 
which  provides  valuable  insights  into  the  nature  and  behavior  of  solu- 
tion trajectories.  In  order  to  keep  this  paper  from  becoming  too  lengthy 
we  restrict  ourselves  herein  to  the  analytic  portion  of  the  analysis 
reserving  the  geometric  analysis  for  a subsequent  paper. 

Our  models  consist  of  two  combat  forces,  an  x-  force  either 
occupying  an  already  supplied  position  or  carrying  enough  supplies  for 
a short  (say  a 30  day)  campaign,  and  a y - force  supplied  by  several 

columns  at  the  levels  s1t...,s  . We  assume  that  the  equations  gover- 

1 P 

ning  the  combat  have  the  form 


* = -aQx  - blV  - b^  bp+1sp  + f(t) 

y - -c1x  - aty  + g(t) 

sk  ■ •Vi*  • Vi\  + rk<t)  ■ 1 ‘ k 1 p • I 


(1) 


In  equations  (1)  the  dots  denote  derivatives  with  respect  to  time.  The 
coefficients  a^  and  a^  represent  losses  in  the  x and  y forces 
not  due  to  direct  fire.  They  may  be  expanded  to  include  such  effects  as 
losses  due  to  desertion,  sickness,  and  accidental  injury  as  well  as  to 
interdiction  by  enemy  fire  power  not  itself  subject  to  attrition.  The 
coefficents  b^  and  c^  are  the  usual  attrition- rate  coefficents 
representing  direct  fire  of  the  combat  forces  upon  each  other.  The  sub- 
system arising  when  b^  “ c^  = 0 , 2 £ k £ p + 1 , has  been  studied  by 
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Bach,  Dolansky,  and  Stubbs  [1]  . 

Our  analysis  includes,  of  course,  complete  results  for  this  subsystem. 


The  coefficients  cj>  2 s j s p + 1 , represent  the  portion 


of  the  direct  fire  of  the  x - force  used  to  interdict  the  supplies  of  the 
y-  force.  We  assume  that  supplying  the  y - force  at  an  adequate  level 
enhances  his  ability  to  inflict  casualties  upon  the  x-  force.  The  coef- 
ficients b2’’",b+l  caPture  this  effect.  The  coefficients  a^  , • • • , 


a^+j  represent  losses  to  the  supply  columns  not  due  to  direct  fire  by 


the  x-  force.  They  may  be  expanded  to  represent  the  effects  of  equip- 
ment breakdown,  adverse  weather  conditions  or  terrain,  as  well  as  losses 
to  fire  power  not  itself  subject  to  attrition.  Finally  the  nonhomogene- 


ous  terms  f(t),  g(t),  r^(t),  1 £ k £ p , represent  reinforcements. 


The  system  (1)  is  to  be  valid  only  for  nonnegative  func- 
tions x,  y,  s^,  1 £ k £ p , which  are  subject  to  certain  stopping  rules 
furnishing  conditions  under  which  one  force  or  the  other  must  disengage 
from  combat.  We  consider  two  types  of  stopping  rules. 


Type  1 stopping  rules:  There  exist  constants  p>0,yi0,ar  >0. 


1 £ j £ p , such  that  the  y - force  must  disengage  if  either 


y S p 


(2) 


or 


y > Qjsj  , 1 * j * P 


(3) 


and  the  x - force  must  disengage  if 


X £ y 


(4) 


Type  2 stopping  rules:  There  exists  cr  > 0 such  that  the  y-  force 
must  disengage  if  either  (2)  holds  or 
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! I 
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and  the  x - force  must  disengage  if  (4)  holds. 

The  condition  (2)  with  0 > 0 is  included  to  avoid  the 
paradoxical  situation  which  would  arise  if  y = 0 and  at  least  one 
8^  > 0 . Then  we  would  have  the  x-  force  suffering  losses  inflicted 
by  supplies  alone. 

The  conditions  (3)  insure  that  the  y - force  is  supplied 
at  an  adequate  level.  They  may  be  regarded  as  the  supply  requirements 
when  the  several  supply  columns  are  furnishing  different  mixes  of  sup- 
plies. The  condition  (5)  should  be  applied  to  the  case  where  each 
supply  column  is  furnishing  approximately  the  same  mix  of  items.  Finally 
the  condition  (4)  is  an  obvious  stopping  rule. 

In  addition  to  the  stopping  rules  we  have  the  initial  con- 
ditions x - Xq,  y * yQ,  Sj  - 3j0  at  time  t = 0 where  all  initial 
values  are  positive  and 

o < * < y0 

yo  s Vjo  ’ 1 S j * p ’ \ (6) 

0 s v < x0 

if  stopping  rules  of  type  1 apply,  and 

2 2 1/2 

o < 8 < y0  s °(sio+'*-  + sPo) 

0 < V < xQ 


if  stopping  rules  of  type  2 apply. 

In  most  of  our  work  we  shall  be  concentrating  on  results 
that  are  primarily  qualitative  in  nature.  This  is  important  because 


t 


either  nonexistent  or  varies  over  such  a large  range  of  values  that  one 
must  conclude  either  that  it  is  unreliable  or  should  be  interpreted  as 
stochastic.  Thus  we  shall  be  interested  in  results  which  depend  only 
upon  the  fact  that  the  coefficents  are  either  positive  or  zero,  but 
which  do  not  depend  upon  their  magnitudes.  Unfortunately,  not  all  ques- 
tions one  might  wish  to  ask  about  solutions  are  independent  of  the  magni- 
tudes of  the  coefficients.  Nevertheless  a remarkable  number  of  properties 
of  solutions  are  qualitative  in  character. 

From  the  system  (1)  we  may  derive  a second  closely  related 
system  in  the  following  way.  Since  we  are  concerned  only  with  positive 
solutions  of  (1)  , let  us  introduce  the  numbers 


and 


a 


max 

lSkSp+1 


Then 


and 


a = min  a . 
ISkSp+l  J 


* * -o*  • Vo  - V! Vi*p  + tu)  ■ 

y £ -c1x  - ay  + g(t)  , 

% * 'Vi*  ' i‘k  + rk(t)  • 1 4 k 1 p • 

* 1 -o'  • Vo  - Vi Vi%  + ■ 

y i -c1x  - ay  + g(t)  , 

*k  * ‘Vi*  • isk  + rk(t)  * 1 S k * p * 


It  follows  that  positive  solutions  of  the  system 


R - - 


5 

• ■ •*ox  ■ Vo  - Vi Visp + f(t> 

Y - -CjX  - ay  + g(E)  (8) 

‘ -W  ■ “k  + rk(t>  • 1 * k s p • 

furnish  upper  bounds  for  positive  solutions  of  the  system  (1)  when 
a * a and  lower  bounds  for  positive  solutions  of  the  system  (1)  when 
a = a . It  is  therefore  of  obvious  importance  bo  carefully  study  pro- 
perties of  positive  solutions  of  the  system  (8)  . 

The  system  (1)  (and,  of  course,  also  the  system  (8))  has  con- 
stant coefficients.  It  follows,  therefore,  that  the  analytic  forms  of 
solutions  are  completely  determined  by  the  eigenvalues  and  eigenvectors 
of  the  coefficient  matrix  A . In  section  2 we  show  how  to  determine 
the  multiplicities  and  the  signs  of  the  eigenvalues  of  this  matrix.  Our 
results  are  complete  and  detailed.  We  also  completely  detail  the  struc- 
ture and  properties  of  the  corresponding  eigenvectors.  Generally  speaking 
all  of  these  results  are  independent  of  the  magnitudes  of  the  elements 
of  the  matrix  A . The  behavior  of  the  largest  eigenvalue  of  A , which 
is  always  simple,  is  an  exception  because  it  can  be  positive,  negative, 
or  zero  depending  upon  the  relative  magnitudes  of  the  elements  of  A . 

Thus  we  can  obtain  both  stable  and  unstable  solutions  of  the  system  (1) 
in  various  cases.  The  system  (8)  also  exhibits  the  same  possibilities. 

Section  3 is  short  and  summarizes  how  solutions  are  explicitly 
constructed  using  the  eigenvalues  and  eigenvectors  of  A . We  do  this 
both  for  the  homogeneous  and  for  the  nonhomogeneous  case.  This  section 
really  consists  of  the  application  of  well-known  results  to  the  present 
problem. 

In  section  4 we  explicitly  construct  the  solution  of  the  system 
(8)  . For  this  special  case  the  eigenvalues  can  be  calculated  as 


functions  of  the  elements  of  A . This  enables  us  to  construct  the 
fundamental  matrix  of  A and  its  inverse.  Rather  simple  analytic 
expressions  are  thus  derived  both  for  the  solution  of  the  homogeneous 
system  and  for  the  nonhomogene ous  system. 

In  this  paper  we  do  not  attempt  to  examine  in  detail  the  role 
played  by  the  stopping  rules.  In  a subsequent  paper  we  will  show  how 
the  stopping  rules  can  be  interpreted  geometrically  as  boundaries  of 
surfaces  which  contain  the  solution  trajectories. 


The  Matrix  of  the  System 


Let  us  set  z “ (x,y ,s^ , • • • ,s  ) and  define  the  matrix 


-0  ‘"I  -b2 


-c„ 


-C, 


0 


1 

0 -a 


•c  .-  0 0 

p+1 


-b 

0 

0 


p+1 


P+ 1 


•J 


(9) 


Then  the  system  (1)  can  be  written  in  the  form 


z = Az  + F(t)  (la) 

where 

F(t)  - <f(t)fg(t),r1(t)>...,r  (t))  . 

Obviously  the  evolution  of  the  solution  of  system  (1)  or  (la)  is 
determined  by  the  properties  of  the  (p  + 2)  x (p  + 2)  matrix  A and 
the  initial  conditions.  In  this  section  we  investigate  the  properties 
of  A . 

The  matrix  A belongs  to  the  class  3^  (see  [2])  which  is  a 


r 
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qualitatively  determined  class  of  matrices  and,  since  symmetrically  placed 
elements  have  the  same  sign,  the  matrix  can  be  symmetrized  (see  [3]).  In 
fact,  for  the  diagonal  matrix 


D+  « diag[l,(b1/c1)1/2,(b2/c2)1/2,...,(bp+1/cp+1)1/2]  , 


we  have  S(A)  = D+*AD+  a matrix  of  the  form 


S (A ) = 


-a. 


.1/2 

(b^)  *(b2C2) 


^1/2 

-(bjCj)  -«1 

-(b2c2)1/2  0 


<ViVi)1/2  0 


■(bp+icp+i) 

0 

0 


1/2 


-a 


p+1 


Thus  the  matrix  A is  similar  to  a symmetric  matrix,  hence  it  has  real 
eigenvalues  and  eigenvectors  and  is  similar  to  a diagonal  matrix. 

We  shall  investigate  more  closely  the  eigenvalues  of  A . To 
this  end  note  that 

p+1  p+1  p+1 

det(XI-A)  = ~f[  (\  + a ) - V b e j]"  (\  + a ) . (10) 

i-0  J-l  J J 1=1 

i/j 

We  now  parallel  the  development  presented  in  the  paper  [4]  (see  also  f5]). 

Suppose  there  are  k distinct  values  among  the  numbers  a^, 

1 s i « p + 1 , which  we  take  to  be  a^,-,',a^  in  increasing  order  and 
suppose  further  that  these  occur  with  multiplicities  m^ , • • • ,m^  res- 
pectively where  obviously  we  must  have  m^+  • • • + * p + 1 . (Note 

that  k * p + 1 is  a possible  value  in  which  case  there  are  no  repeti- 
tions among  the  numbers  a , 1 t i i p + 1 ,)  From  the  righthand  side 

1 k ■ -1 

of  (10)  we  can  factor  out  | f (X.  + a^)  , hence  -a^  is  an 


. . 
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eigenvalue  of  A of  multiplicity  m^  - 1 . Now  divide  det(xI-A)  by 


m 

- si 


1 1 (X  + a ) to  obtain 
i=l 


v»he  re 


det(XI  - A)..  „ + _ 

k m X 0 X 

TTu+V 

i=i 


k d 

" dft)  + *o  ■ 


k d . 

dM  ■ 1 • £,r?r 

1=1  J 


and  d is  the  sum  of  the  m.  values  of  b.c.  associated  with  a.  . 

j J 1 1 J 

Clearly,  the  eigenvalues  of  A not  found  among  the  a,  , 1 s j i k , 
are  going  to  satisfy 


d(X)  = -a. 


(ID 


Let  us  plot  d (x)  against  \ . Clearly,  d(X)  is  asymptotic  to 
each  of  the  vertical  lines  X = a j > 1 s j s k . Moreover,  we  have 


d'(X>  - 1 + Yj  L- ? > 0 , 

j=l  (X  + a ) 

hence  d is  an  increasing  function  of  X wherever  it  is  defined. 

Also  d(0)  <0  if  all  a^  > 0 and  d(0)  is  not  defined  if  a^  = 0 , 
d -*  oo  as  X **  00  and  d -»  -«  as  X -»  -®  . It  follows  that  the  graph 
is  as  shown  below. 


\ -'.m 


$ 


a iafS  dkUaa 


0 the  vertical  axis  is  a vertical  asymptote.) 


We  see  that  the  k + 1 roots  of  (11)  , denoted  by  p, 


the  inequalities 


It  follows  that  the  p + 2 eigenvalues  of  A are  divided  into  two 


groups  as  follows 


I.  There  are  p + 1 - k eigenvalues  consisting  of  m 


which  are 


II.  There  are  k + 1 eigenvalues  pi 


The  first  group  can  be  empty  (the  case  k = p + 1 ) but  there  are 


always  at  least  two  elements  in  the  second  group.  In  the  case  k = 1 


the  sequence  of  eigenvalues  has  the  form 


with 


occurring  p times 


10 


"■ *"  " ■ 


_ I — 


It  is  important  to  observe  that  the  largest  eigenvalue  ^ can 

be  positive,  negative,  or  zero  (the  case  > 0 is  illustrated  in 

figure  1).  Obviously,  if  a^  = 0 or  a^  = 0 then  > 0 , but  if 

a^  < 0 and  a^  > 0 all  cases  can  occur.  Moreover,  we  can  formulate 

the  conditions  for  X = 0 to  satisfy  (11)  . In  fact,  d(0)  = 
k d 

- } T-*  , so  p,  = 0 is  an  eigenvalue  if 

j=laj 

k d. 

3°  = 

j=i  i 


(13a) 


< 0 if 


and  ^.  > 0 if 


•o>  if*  • 


j=i  j 

k d 


a0<  E J 


j-laj 


(13b) 


(13c) 


In  the  case  where  >j,  < 0 the  matrix  A is  a stable  matrix,  otherwise 
A is  not  stable.  We  shall  see  subsequently  that  the  behavior  of  solutions 
of  the  system  (1)  are  not  qualitatively  different  in  the  stable  case 
from  their  behavior  in  the  nonstable  case. 

It  is  instructive  to  examine  two  specific  examples  of  the  class 
of  matrices  A . 

Example  1.  Consider  the  case  p = 1 so  that  (Sq>0) 


1 

j 

A 85 

-a  ^ -b  -b 

I 

) 

J 

0 1 2 

-c,  -a,  0 

1 

\ 

1 1 

-c_  0 -a 

t 

! 

P 

fSJ 

} 

We  have  det(Xl-A)  “ (X  + aQ)  (X  + ) (X  + a? ) - b^U  + a^  - b^U  + a^  . 

Thus  two  basic  cases  can  arise.  First,  if  a^  = ■ a^  , we  have 


likiitieiittJifelia 


wff&r:  :::± 


^ 


n 


2 blCl+b2C2 

det(XI-A)  - (X+ij)  [X+aQ j^jT| 1 • Tt  follows  that  - a ] is  a 


simple  eigenvalue  of  A whose  three  eigenvalues  satisfy 


M.2  < <M-1  • 


For  a^  = 0 , p,^  > 0 . Otherwise  A is  a stable  matrix  if 


‘o’ 


Vl+b2C2 


b c +b  c 

and  A is  not  stable  if  a_  £ = . Next  if  a f a and 

u a-^  L L 


min(a^,a2)  < 0 , we  have 


-p,  < -tnax(a1,a2)  < p*2  < -minCa^a^  < p^  . 


Now  A is  a stable  matrix  if 


_ Vl  ^ b2C2 

> + 


0 a^  a2 


and  A is  not  stable  if 


* Vl  * Vj 

a„  s — — + 


0 al  a2 


Note  that  multiple  eigenvalues  cannot  occur  when  p = 1 . 


Example  2. 
We  have 


Consider  next  the  case  p “ 2 (again  assuming  a^  > 0 ) 


‘o  -bI  "b2 


-b 


-c. 


- c. 


1 

0 -a. 


0 

0 

-a„ 


! 


so  that 


12 


det(\I-A)  « (\  + a0)(\+ai)(X+a2)(^  + a3)  - b c^X  + a2  ) (X  + ) 

“ b2c2^  + aj)^  + a3)  " b3C3^  + a1)(^  +a2>  • 


There  are  now  essentially  three  cases. 


Case  1.  a^  “ a2  * a^  = a^  > 0 . Then 

i be  + be. + b.c„ 

det(\I  -A)  - (X+a1)3[X+aQ  - 1 1 ^ +2fl2 2_J] 

and  the  eigenvalues  of  A satisfy 


^2  < "al  < ^1 


with  -a^  occurring  twice  as  an  eigenvalue.  A is  stable 


if 


ao> 


b, c,  + b c„  + b„c 
1 1 2 2 3 3 


and  nonstable  if 


aos 


blCl  + b2C2+b3C3 


i 


Case  2.  Two  of  the  numbers  a^  a2  , a are  equal  and  the  third  different, 


say,  a3  - a3  = a2 ’a2  = a\  * Now 

2 bici+b9C9  b9C-j 

det(XI-A)  = (X+i  )(X  + i,r[X+an  - - } - T^ 

1 2 0 X+a2  X + a^ 


and  the  eigenvalues  of  A satisfy 


1 


$ 


< -a3  <*2  < -ix 


with  -a2  a simple  eigenvalue.  This  time  A is  stable  if  and  only  if 


blc1+b3c  b9c9 
a„  > . — — + -i-=- 

0 a„  a. 


{ 

1 
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Case  3.  The  numbers  a^,  are  distinct  with,  say,  a?  < < a.( 

hence 

Vl  b2°2  b3c3 

det(XI  - A)  - a + «1)tt  + a2)a  + «3HX+*0-^I-  - 5^  - ^1 
so  the  eigenvalues  of  A satisfy 


^4  < a3  < ^3  < 81  < ^2  < a2  < *1  ' 

blCl  b2C2  b3C3 

Now  A is  stable  if  and  only  if  & > + + 

1 0 aj  &2  a j 


Let  us  turn  next  to  the  eigenvectors  of  A . We  must  consider  the 
equation  (XI-A)^  = 0 where  u **  (ui>’’->un)T  and  X is  an  eigenvector 
of  A . This  system  is 


&1 

(\  + a )u.  + ) b u 

01  j«l  j 


J+l 


- 0 , 


cjUl  + (X  + aj)uj+1  = 0 . 


(14) 

(15) 


Suppose  first  that  X “ , 1 £ j s k + 1 , then  X + f1  0 foi 

1 & j £ p + 1 , and  from  (15)  we  obtain 


l-l 

u<  TZ ui  * 2 s i s p + 

1 Vi-i1 


Choosing  Uj  * 1 , we  obtain  the  eigenvectors 


c c 

“j  ■ 


^ . + a 


(16) 


J pH-1 


1 s j s k + 1 . Now  c^>0  , lSiSp+1  , hence  the  signs  of  the 
elements  of  the  eigenvectors  depend  only  upon  the  denominators  of  the 
components.  For  J **  1 we  have  p. ^ + a^  > 0 for  1 s i S p + 1 , 
hence  every  component  is  negative  except  for  the  first  component.  Thus 
"u^  has  the  sign  pattern  . On  the  other  hand  for 


,.T* 


J 
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k+1 


J “ k + 1 we  have  p^+j  + a^  < 0 for  IsiSp+l  , so  that  u 

has  the  sign  pattern  (+,+,+,••♦,+)  . Finally,  consider  , 

1 < j < k + 1 . We  have  + a^  < 0 for  + •••  + ^ values  of 

a^  and  + a^  > 0 for  m^  + values  of  a^,l£i£p+l. 

Thus  each  component  of  u^  for  which  a^  has  one  of  the  values 

^ will  be  positive  and  each  component  for  which  a^  has  one 

of  the  values  a.,*'*,3  will  be  negative. 

J k 

Next  suppose  \ = -a^  for  some  1 £ i £ k . Then  m^  equations 
of  the  system  (15)  reduce  to  the  form  c^u^  = 0 » hence  = 0 . 

The  remaining  p + 1 - m^  equations  from  this  subsystem  thus  reduce  to 
(-a^+a^)u^+i  = 0 . Since  (-a^+a^)  f 0 , it  follows  that  the  corres- 
ponding Uj  = 0 . Finally  the  m^  components  which  need  not  be  zero 
must  satisfy  (14)  , which  becomes 

p+1 


ft 


biVi  ’ 0 


(14.) 


because  u^  ■ 0 . There  are  m^  nonzero  terms  on  the  lefthand  side 
of  (14a)  , and  we  can  clearly  obtain  m^  - 1 linearly  independent 
solutions,  hence  m - 1 linearly  independent  eigenvectors.  In  fact, 
it  is  easy  to  see  each  of  these  eigenvectors  can  be  chosen  so  as  to 
have  exactly  two  nonzero  components,  one  positive  and  one  negative. 

In  order  to  see  exactly  how  the  eigenvectors  appear  in  the  various 
possible  cases  let  us  find  them  in  the  special  case  p = 2. 


Example  2 (continuation).  Consider  first  the  eigenvectors  for  case  1, 

a. 


1 2 3 

eigenvectors 


il>0 


Corresponding  to  p.^  and  p,^  we  have  the 


(1,  - 


fti+al 


+ a 


) , i - 1.2  , 


while  corresponding  to  the  double  eigenvalue  -a^  we  have  the  eigenvectors 


u3  = (0,b2,-br0)  , 


U4  = (°*b3»0*'b1^  • 


Turning  next  to  the  case  2,  we  have  the  eigenvectors  u^,  u2 , u^ 


given  by 


C1  C2  C3 

u.  = (1, ■£=-  , -■  ) , i = 1,2,3 

i M-  i + a ^ +a  m,  +a 


and  the  eigenvector 


u4  = (0,b;j,0,-b1)  . 


Finally,  in  case  3 all  eigenvectors  u^,**’,u4  have  the  form 


- (1, 


^i  + al’  ’^i  + V ’‘ii  + a3 


) , 1 s i s 4 . 


J 

1 


We  wish  to  point  out  that  one  of  the  consequences  of  our  analysis 
to  this  point  is  that  multiple  eigenvalues  can  occur  only  in  the  event 
that  there  are  two  or  more  supply  columns.  This  means  that  the  matrix 
A must  be  at  least  of  order  4. 

It  is  of  considerable  interest  to  understand  how  the  eigenvalues 
of  A vary  when  the  elements  of  A vary.  To  this  end  consider  first 
the  eleirent  a^  . It  is  clear  from  figure  1 that  increasing  aQ  will 
decrease  each  of  the  eigenvalues  ^ , 1 s ] s k + 1 . Thus  we  can 
assert  that 


<0,  lsjsk+1. 


On  the  other  hand,  it  is  also  clear  that  the  eigenvalues  -a^  , isiik  , 
do  not  depend  upon  an  , i.e., 


^ ft 


^■fll 


m irfi 


MRR0VM 


d(-ai) 

da„ 


0 , 1 £ i s k 
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(18) 


Thus  we  have  complete  information  regarding  the  dependence  of  the  eigen- 
values upon  the  element  a^  . 


3.  Construction  of  Solutions 

Now  it  is  well  known  (see  [6],  [7])  that  the  solution  of  the  homo- 
geneous system  (la)  can  be  written  in  the  form 


z(t)  = eWz(0)  (19) 

T 

where  z(0)  ■ * * " ,8p0^  vec^or  °i  initial  values  at 

time  t * 0 . The  matrix  exp(tA)  is  called  a fundamental  matrix  of 
the  system.  Let  us  determine  the  form  of  this  matrix. 

To  this  end  let  U be  the  matrix  whose  column  vectors  are  the 
eigenvectors  of  A . Then 


AU  - UJ 

where  J is  the  Jordan  matrix  of  A . Since  A is  similar  to  a 
symmetric  matrix,  J is  a diagonal  matrix  with  diagonal  elements  the 
eigenvalues  , and  -a^,***,-a^  , the  latter  occurring  with 

appropriate  multiplicities.  Thus,  denoting  the  eigenvalues  of  A by 
V-.X*  . have 

J - diagUj,  • • • , \nl  . 


Now  we  have 
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tl  \ t 

e = diag[ e 1 , ' • * ,e  n ] . 


More  explicitly,  then 


tA  “l*  ^k+l*"  _alt  'V.  -1 

e = U diag[e  ,e  .•••,e  ]U 


It  follows,  therefore,  that  (19)  appears  in  the  form 


^l1-  -1 

z(t)  = U diag[ e .---.e  K )U  z(0)  . 


(19') 


Let  us  examine  this  formula  more  closely.  First  consider  the 
^’lt  '*kt 

matrix  U diag[e  , •••,e  ] . Denote  this  matrix  by  U(t)  . Multi- 

plication of  any  matrix  on  the  right  by  a diagonal  matrix  multiplies  the 
columns  of  the  matrix  by  the  corresponding  diagonal  elements.  Thus  the 
first  k + 1 column  vectors  of  U(t)  are  just  the  vectors 


u,(t)  = (e  J , - j.  . , • ‘ • , 


J\  ~ w | i _ f y _l_- 

“j+*i  “j+Vi 

1 s j s k + 1 . The  remaining  columns  of  U have  the  form 


(20) 


-v 


a . t 
l 


V0  " (“ki®  ’ * * ’ ,akp+2e  > 


(21) 


where  exactly  two  of  the  numbers  , 1 £ j s p + 2 are  different 

from  zero  and  each  a^  is  one  of  the  numbers  a^,‘‘*,a^  . 

Next  we  observe  that  U *z(0)  is  a vector.  Let  us  denote  it  by 
w . Then  the  formula  (19)  can  be  rewritten  as 


z(t)  * U(t)w 


(22) 


Let  us  turn  next  to  the  solution  of  the  nonhomogeneous  system 


z “ Az  + F(t) 


(23) 


where  f(t)  ■ (f(t) ,g(t),r^(t) , • • . ,r^(t))  . As  mentioned  above  we  may 
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Chink  of  the  vector  F as  representing  reinforcements.  In  the  case  of 
the  x and  y forces  this  might  be  additional  personnel,  while  for  the 
supply  columns  it  might  consist  of  vehicles  or  other  equipment  replace- 
ments. The  general  formula  for  solving  the  system  (23)  is 

z(t)  - e^O)  + Jo  C e(t's)AF(s)ds  . (24) 

In  view  of  the  form  which  we  have  presented  above  for  the  fundamental 
matrix  of  the  system,  we  can  rewrite  (24)  as 

z(t)  = U(t)w  + f U(t  - s)u’1F(s)ds  . (25) 

0 

It  is  evident  from  the  formulas  (22)  and  (25)  that  we  must 
investigate  the  form  of  the  matrix  U * . We  shall  not  attempt  to  do 
this  for  the  most  general  case. 


4 . Upper  and  Lower  Solutions 

As  was  pointed  out  in  section  1 the  system  (8)  is  of  particular 
importance  because  its  solutions  for  appropriate  values  of  a can  be 
used  to  obtain  both  upper  and  lower  bounds  for  all  positive  solutions 
of  the  system  (1)  . In  this  section  we  shall  explicitly  solve  the 
system  (8)  . 

Thus  consider  the  matrix 


! 

j 

A 8 

Ao 

-aQ  -b1  . . . -b^ 

1 

-Cj  -a  ...  0 

• • • • 

• • • • 

-c  . - 0 ...  -a 

l’11  J 

The  characteristic  polynomial  of  Aq  Is 


i 


tti'i  tii 


det(XI  - Aq)  = (\+a)r[x  + (*Q  +a)X  +aaQ  - 6 ] , 

P+1  / 

where  6 “ Jj1bjCj^  * Therefore  the  eigenvalues  of  AQ  are 


X = -a  , p - times  , 


X=  +][(a()-a)2  +462lV/2 


X = - 


ao  + a 1 


2 |f(a0  - a)2 +4<52]  1//2  . 


Let  us  set 


|/(a0-a)2+4fl2  , 

ao  + a 


so  that  the  eigenvalues  of  Aq  becone 


X - -a , p - times  , 


X ° a + e 


X = ct  - e . 


Now  the  fundamental  matrix  corresponding  to  A Is 


1 0 ...  0 

C1  b2  • ’ * bp+l' 

o - e +a  -b1  . . . 0 


a + E + a cr-E  + a 


It  Is  easy  to  show  that  det  U * (- 1 )P2bP" *e  and  one  obtains 


■ .. 


where  we  have  defined 


Next  let  ua  define 


and  compute  the  components  of  the  vector  w.  After  some  calculations 


one  finds 


With  the  vector  w at  hand  we  can  write  out  the  solution  to 


the  homogeneous  version  of  the  system  (8).  It  is 


x ( t ) - e°  {xQ  cosh  et  + -[ (o  + a)xQ  - DQ] sinh  et}  , 
'lD0,  -t  Cl.  to 


c.D_  . c.  (ct  + a)D~  . c,Drt 

,.v  , 1 0,  -at  lr  . 0,  ot  , , 1 0 at 

y(t)  ° (yQ 2~ )e  - — lxo+ 2 *e  et  + — 2~e  cosh  et> 

6 6 fi 


, Ck+1D0.  -at  Ck+lr 

sk(t)  = (sk,0 ^2  e ' ~rfx0  + 


(a  + a)D 


-]eatsinh  et 


Ck+1D0  at  . _ . 

+ 5 — e cosh  et  , 1 £ k £ p . 


Now  observe  that  the  stable  case  occurs  if  and  only  if 


aa  - 6 >0 

0 


But  we  also  note  that  in  any  event  a ^ 0,  a < 0,  and  e > 0 . Conse- 
quently the  analytic  form  of  the  solution  (28)  is  the  same  whether  or 
not  Aq  is  stable.  This  is  attributable  to  the  fact  that  only  the  sign 
of  the  eigenvalue  \ = a + e is  affected  by  the  change  from  the  stable 
to  the  unstable  case.  However,  we  also  wish  to  call  special  attention 

to  the  situation  occurring  when  a = 0 . Then  the  solution  will  always 

2 

be  unstable  because  (29)  reduces  to  -6  >0  which  can  never  be 

satisfied. 


Let  us  next  write  out  the  solution  to  the  nonhomogeneous  problem. 
To  this  end  it  is  sufficient  to  write  out  the  integral  terms  appearing 
in  equation  (25)  because  these  terms  added  to  the  corresponding  terms 
in  (28)  yield  the  entire  solution. 


Set  w = U F(s)  . Then  we  have 


”2  + 
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denote  the  solution  of  the 


nonhomogene ous  system  satisfying  zero  initial  conditions.  We  then  obtain 


{f(s)cosh  e(t  - s)  + [ 


Again  because  a a 0 , a < 0 , and  e > 0 this  portion  of  the 


solution  has  the  same  form  in  both  the  stable  and  the  nonstable  case 


The  formulas  (28)  and  (31)  provide  the  basic  tools  for  any 


additional  analytic  analysis  of  the  system  (8).  On  the  other  hand,  as 


was  pointed  out  in  the  introduction,  the  solutions  we  have  obtained 


provide  upper  and  lower  bounds  for  all  solutions  of  the  system  (1) 


These  bounds  will  hold  as  long  as  the  bounding  solutions  remain  positive 


and  do  not  violate  the  stopping  rules.  Since  the  bounding  solutions 


have  the  same  Initial  values  as  the  solutions  we  wish  to  bound,  either 


the  upper  bound  or  the  lower  bound  may  be  expected  to  remain  in  force 


. dLUrW*  KJJJir.  ' T:.'  ■ ■•.T" 


until  the  solution  violates  one  c che  stopping  rules.  In  view  of  the 
fact  that  both  bounding  solutions  have  the  same  mathematical  form  we 
expect  that  they  will  furnish  reasonably  close  estimates  of  the  actual 
solution.  This  is,  of  course,  only  a conjecture.  We  expect  to  examine 
more  closely  the  behavior  of  the  solutions  of  our  model  in  a subsequent 


paper. 


In  view  of  the  fact  that  our  aim  in  this  paper  was  to  initiate 


the  study  of  models  whose  dependence  upon  the  model  parameters  is 
primarily  qualitative  it  is  gratifying  to  observe  that  the  upper  and 
lower  solutions  both  have  the  same  structure  independently  of  the  magni- 
tudes of  the  elements  of  A as  long  as  all  of  the  numbers  a,,  , 0 £ i £ p 
+1,  b^  , 1 £ i £ p + 1 , are  positive.  This  is  about  the 

strongest  qualitative  result  we  could  hope  to  obtain.  Moreover,  the 
form  of  the  solution  in  the  case  where  a^^O,  0 £ i £ p + 1 is 

also  easily  deduced  from  the  above  results. 


" 'i 


. . - 
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